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Abstract 

In this paper we study the scahng behavior of the interface fluctuations (roughness) for a discrete 
model with conservative noise on complex networks. Conservative noise is a noise which has no 
external flux of deposition on the surface and the whole process is due to the diffusion. It was 
found that in Euclidean lattices the roughness of the steady state Wg does not depend on the 
system size. Here, we find that for Scale-Free networks of N nodes, characterized by a degree 
distribution P{k) ~ /c"^, Ws is independent of N for any A. This behavior is very different than 
the one found by Pastore y Piontti et. al [Phys. Rev. E 76, 046117 (2007)] for a discrete model 
with non-conservative noise, that implies an external flux, where Wg ~ InA^ for A < 3, and was 
explained by non-linear terms in the analytical evolution equation for the interface [La Rocca et. 
al, Phys. Rev. E 77, 046120 (2008)]. In this work we show that in this processes with conservative 
noise the non- linear terms are not relevant to describe the scaling behavior of Wg. 

PACS numbers: 89.75.Hc 68.35. Ct 05.10.Gg 05.45.Xt 
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I. INTRODUCTION 



It is known that many physical and dynamical processes employ complex networks as the 
underlying a substrate. For this reason many studies on complex networks are focused not 
only in their topology but also in the dynamic processes that run over them. Some examples 
of these kind of dynamical processes on complex networks are cascading failures traffic 
flow , epidemic spreadings [4] and synchronization |5|, |6|]. In particular, synchronization 
problems are very important in the dynamics and fluctuations in task completion landscapes 

FT 

in causally constrained queuing networks [7\, in supply-chain networks based on electronic 
transactions 8^, brain networks js], and networks of coupled populations in the synchroniza- 
tion of epidemic outbreaks [lO|. For example, in the problem of the load balance on parallel 
processors the load is distributed between the processors. If the system is not synchronized, 
few processors have low load and they will have to wait for the most loaded processors to 
finish the task. The nodes (processors) of the system have to synchronize with their neigh- 
bors to ensure causality on the dynamics. The computational time will be given by the most 
loaded processors, thus synchronizing the system is equivalent to reduce or optimize the com- 
putational time. Synchronization problems deal with the optimization of the fiuctuations 
of some scalar field h (load in processing) in the system that will be optimal synchronized 
minimizing those fiuctuations. To analyze synchronization problems is customary to study 
the height fiuctuations of a non-equilibrium surface growth. If the scalar field on the nodes 
represents the interface height at each node, its fiuctuations are characterized by the average 
roughness W{t) of the interface at time t, given by W = W{t) = \l/NJ2iLi{hi - (/i))^}^''^ , 
where hi = hi{t) is the height of node i at time t, (h) is the mean value on the network, N 
is the system size, and {.} denotes an average over configurations. Pastore y Piontti et. al 



ll| studied this mapping in Scale-Free (SF) networks 12|] of broadness A and size using 



a surface relaxation growth model (SRM) [13] with non- conservative noise and found that 
for A < 3 the saturation roughness Ws scales as Wg ~ IniV. Later, the evolution equa- 
tion for the interface in this model was derived analytically [l^ for any complex networks. 
The derived evolution equation has non-linear terms as a consequence of the heterogeneity 
of the network that together with the non- conservative noise are necessary to explain the 
Ws ~ InA^ behavior for A < 3. However, there exist many physical processes where the 
noise is conservative and cannot be modeled as a fiux deposition on a surface. In models 
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without external flux, where particles are moved by diffusion, the total volume of the system 
remains unchanged. Examples of this kind of process are thermal fluctuations, diffusion by 
an external agent such as an electric fleld, load balance of parallel processors where the total 
load in the system is constant over a certain time interval. For the last example, the only 
flux is due to diffusion of the load from a processor to another. Though not extensively. 



conservative noise has already been studied in Euclidean lattices 
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16| and it was found 



that Ws does not depend on the system size L. The evolution equation of this process in 



Euclidean lattices can be well represented by an Edwards- Wilkinson (EW) process [17| with 
conservative noise. In this paper we study this model in SF networks by simulations of the 
discrete model (Section [IT]) and derive analytically its evolution equation for any complex 
networks (Section IIIip . Those networks represents better the heterogeneity in the contacts 
in real systems, like the Internet, the WWW, networks of routers, etc.. We applied the 
mean-fleld approximation to the evolution equation and show that the scaling behavior of 
Ws with N (Section HVT) is only due to flnite size effects. To our knowledge this class of 
model was never studied before in complex networks. 

II. SIMULATIONS OF THE DISCRETE MODEL 

In this model, at each time step a node i is chosen with probability If we denote by 
Vi the nearest-neighbor nodes of i , then (1) if hi < hj, Vj G f j ^ the scalar flelds remains 
unchanged, else (2) if hj < hn, Vri ^ j & Vi ^ hj = hj + 1 and hi = hi — 1. In that 
way the total height of the interface is conserved and we have that the average height is 
constant. We measure the roughness W for SF networks, characterized by a power law tail 
in the degree distribution P{k) ~ A;^^, where kmax > k > k^in is the degree of a node, k^ax 
is the maximum degree, kmin is the minimum degree and A measures the broadness of the 
distribution [l2]. To build the SF network we use the MoUoy Reed (MR) [l8| algorithm or 
configurational model. 

In Fig. [1] we plot W"^ as a function of the time t for different system sizes and in Fig. [2] 
the steady state as a function of A^, for (a) A = 3.5 and (b) A = 2.5. We can see that 
increases with but, as we will show later, this dependence in the system size is only due 



to 



inite size effects introduced by the correlated nature (dissortative) of the MR algorithm 



19| . For all the results we use kmin = 2 in order to ensure that the network is fully connected 
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20| . and assume that the initial configuration of {hi} is randomly distributed in the interval 
[—0.5,0.5]. Then, we have that (h) = 0. 



III. DERIVATION OF THE STOCHASTIC CONTINUUM EQUATION 



Next we derive the analytical evolution equation for the scalar field hi for every node i in 
the conservative model in random graphs. The procedure chosen here is based on a coarse- 
grained (CG) version of the discrete Langevin equations obtained from a Kramers-Moyal 



expansion of the master equation 
of the height in any growth model is given by 

dhi 1 



2lN23|. The discrete Langevin equation for the evolution 



22 



23| 



dt 



T 



(1) 



where takes into account the deterministic growth rules that produces the evolution of 
the scalar field hi on node i, t = N6t is the mean time of attempts to change the scalar 



fields of the interface, and rji is a noise with zero mean and covariance given by 22|, |23| 



T 



(2) 



More explicitly, and Kfj are the two first moments of the transition rate and they are 
given by 

N 

(3) 



Kl — ^ Aij [Pij — Pji] , 



N 



K 



6ij ^ Ain{Pin + Pni){^ 



5. 



I] J 5 



(4) 



n=l 



where {^ij} is the adjacency matrix {Aij = 1 if i and j are connected and zero otherwise) 
and Pij is the rule that represents the growth contribution to node i by relaxation from its 
neighbor j. In our model the network is undirected, then Aij = Aji. As the rules for this 
model are very complex if we allow degenerate scalar fields, we simplify the problem taking 
random initial conditions [See discrete rules on SecIU]. Thus, 

p^^ = Q(^h, - hi) n [1 - Q{K - hn)] , 

where B is the Heaviside function given by 0(x) = 1 if a; > and zero otherwise, with 
X = ht — hs = Ah. Without lost of generality, we take r = 1. 
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In the CG version Ah — )■ 0; thus after expanding an analytical representation of 0(x) in 
Taylor series around a; = to first order in x, we obtain 

N N 

i=i i=i 

N 

^jnihn ~ hi) 



C C ^ 



0((A/.)^) , (5) 



where Cq and Ci are the first two coefficients of the expansion of 6(x) and Q{ki) = (1 — Co)'^'^'^ 
is the weight on the link ij introduced by the dynamic process. 

Notice that the network is undirected and the noise is conservative, thus the average 
noise correlation [see Eq. is {rii{t)rij{t')) = 0, where ( ) represents average over all the 
nodes of the network. Notice that in Eq. (jS]) the non-linear terms are disregarded. As we 
will show below, for this conservative noise model these terms are not necessary to explain 
the scaling behavior of Wg with A^. 

We numerically integrate our evolution equation Eq. ([1]) in SF networks using the 
Euler method with a representation of the Heaviside function given by Q{x) = (1 + 
tanh[f/(x + z)])/2, where U is the width and z = 1/2 23|]. With this representation, 
Co = (l + tanh[[//2])/2 and Ci = {1 -tanh^[U/2]) U/2. We assume that the initial configura- 
tion of {hi} is randomly distributed in the interval [—0.5, 0.5] and for the conservative noise 
we used the algorithm described in {2^: at each time step, for every node in the network 
and for any of its nearest neighbor we add a random number in the interval [-0.5,0.5] and 
remove this amount to one of the nearest neighbor nodes. 

In Fig. |3]we plot as a function of t from the integration of Eq. ([T]) for (a) A = 3.5 
and (b) A = 2.5, and different values of N with kmin = 2. For the time step integration 
we chose At <^ ^/Knax according to Ref In Fig. IHwe plot the steady state as a 
function of for (a) A = 3.5 and (b) A = 2.5. We can see that depends weakly on A^, 
but as shown below this size dependence is due to finite size effects introduced by the MR 
construction. Next we derive the mean field approximation in order to explain the nature 
of the corrections to the scaling. 
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IV. MEAN FIELD APPROXIMATION FOR THE EVOLUTION EQUATION 



We apply a mean field (MF) approximation to the linear terms of Eq. ([5]). In this ap- 
proximation we consider 1 ^ /c^m ^ k^ax- Taking Cij = AijQ{kj) and Tjj„ = AijAjn^{kj), 
then 



Kl = Co [Q - hn{k,)] + c,Q 



+ c\VL{ki)k, 



+ 



where 



:i - Co) 

Ti -- 



^ A-h- 
j=i '^^ 

^ Af AT T h 

E ^ 



(6) 



j = l n=l,nj^i 



rp 



1=1 Z^n=l 



(7) 



Disregarding the fluctuations, we take Y.f=i^ijhj/ki ~ (/i), Y.f=iCijhj/Ci ~ (/i), and 
J2f=iJ2n=i,njti'^ijnhn/Ti ^ (/i). From Eq. (I7j), we can approximate d by Ci{ki) ^ 
ki J^^°^ P{k\ki) f2(/c) d/c jsl, where P(/c|fcj) is the conditional probability that a node with 
degree fcj is connected to another with degree k. For uncorrelated networks P{k\ki) = 
kP{k)/{k) Q, then Ci{ki) ^ h ki/{k) with Ji = jj^^^^ P{k) k n{k) dk . Making the same 
assumption for Tj, we obtain Ti{ki) ^ I2 h/ (k) with Is = J^27 ^i^) k {k - 1) ^{k) dk. 
Then, the linearized evolution equation for the heights in the MF approximation can be 
written as 



dhj 
dt 



Fi{ki) + Ui{h) {{h)-hi)+7]. 



(8) 



where Fiiki) = c^ki [Ii/ {k) — Q{ki)] represents a local driving force, Uii^ki) = cikiip + 

is a local superficial tension-like coefficient with h = + /2Co/(l — cq)]/ {k). This mean field 

approximation reveals the network topology dependence through P{k). 

Taking the average over the network in Eq. (|8]), d{h)/dt = l/A^I]^i Fi = 0, then (h) is 



constant in time. The solution of Eq. (|8]) 2l| is given by 



hiit) 



Fi + Uiih) 



-''^^'-'^r]i{s) ds . 



(9) 



Using Eq. (Q and the fact that in our model with the initial conditions we use (h) = 0, 
we find the two-point correlation function 

mi)h,{h)} = (^) (^^^ (1 - e--*)(l - e-^') 
Jo Jo 

For t > max {l/z/j}, we can write Wg as 

^ {<.?>} 4i: f .lEf, (10) 

1=1 ' 1=1 ' 

where Kf^ [See Eq. (j4])] is given by 



TV r 



i=i 



For SF networks it can be shown that /i,/2 ~ const. + /cmax exp(—fcmax const.), where 
kmax ~ A^i/(^-i) for MR networks; thus we can considerer the quantities Ji and I2 as 
independent of A^. 

From Eq. (fTOl) and using the expressions for Fj, i/j and i^^fj, we have 

'^^-(|)^B'^^g(/-ft))^ + ^B^E/.W. (11) 

where 

1 ± 



and -8 = 1/(1 + Co/2/ (1 — Co)/i). Taking the continuum limit we find another expression for 
Eq. (dH) as 

W^, = [^) B' p(k)if4k)y dk+^B p{k)U{k)dk. 

The function f±{k) has a crossover a.t k = k* , where k* is the crossover degree between 
the two different behaviors, then 

1) for k<k* ^ f±{k) ^ ± {k) n{k)/h/2, and 

2) for k> k* ^ f±{k) ^ 1. 

As k* is the crossover between two different behaviors of f±{k), and the numerator of 
the function diverges faster than the denominator, we have 1 ~ (k) Q{k*)/Ii thus k* ~ 
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ln(Ji/(fc))/ln(l-Co). Then, 



\CiJ Ci Jk* 



p{k)dk + 



+ 




Even thought k* depends on k^ax, it can be demonstrated that the two last integrals 
depends weakly on and can be considerer as constant. Then, introducing the corrections 
due to finite size effects trough k^ax in {k), we obtain 



where Ai, A2 and As do not depend on kmax- 

In Fig. |2] and Fig. H] the dashed lines represent the fitting of the curves with Eq. ( IT2|) 
considering finite-size effects introduced by the MR construction. We can see that this 
equation represents very well the finite size effects of this model. This means that even 
though the networks is heterogeneous, the non-linear terms are not necessary to explain 
the N independence of Wg when a conservative noise is used. Notice that even when our 
network is correlated in the degree, the expression for found describe very well the 
scaling behavior with N as shown in the insets of Fig. [2] and HI This model suggest a useful 
load balance algorithm suitable for processors synchronization in parallel computation. Our 
results show that the algorithm could be useful when one want to increase the number of 
processors and its general behavior is well represented by a simple mean field equation. 

V. SUMMARY 

In summary, in this paper we study a conservative model in SF networks and find that 
the roughness of the steady state is a constant and its dependence on for any A it is 
only due to finite size effects. We derive analytically the evolution equation for the model, 
and retain only linear terms because they are enough to explain the scaling behavior of Ws- 
Finally, we apply the mean field approximation to the equation and we calculate explicitly 
the corrections to scaling of Ws- This approximation describe very well the behavior of the 
model and shows clearly that the corrections are due to finite size effects. 




(12) 
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256 (o), 512 (A), 1024 (y), 2048 (+), 3072 (*) and 4096 (X) and b) A = 2.5 for = 64 (o), 
128 (□), 256 (o), 512 (A), 768 (y), 1024 (+) and 1280 (*). Each curve was obtained with 10.000 
reahzations. 
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FIG. 2: Wg as a function of N for a) A = 3.5 and b) A = 2.5 in symbols for the same system sizes of 
the Fig [TJ The dashed Unes represent the fitting with Eq. (jl2p , obtained in the MF approximation 
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FIG. 3: W'^ as a function of t from the integration of the evokition equation: a) A = 3.5 for = 64 
(o), 128 (□), 256 (o), 512 (A), 1024 (y), 1536 (+) and 2048 (★) and b) A = 2.5 for iV = 64 (o), 
128 (□), 256 (o), 512 (A), 768 (y), 1024 (+) and 1280 (★). For all the integrations we use U = 0.5 
and typically 1.000 realizations of networks. 
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FIG. 4: W'l as a function of iV for a) A = 3.5 and b) A = 2.5 in symbols for the same system 
sizes of the Fig [3l The dashed hnes represent the fitting with Eq. (112p , obtained by considering 
the finite-size effects introduced by the MR construction. 
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